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Abstract 

The  general  aim  of  this  project  is  to  get  an  improved  understanding  of  the  interac¬ 
tion  between  wall-generated  turbulence  and  compliant  surface  coatings  using  analysis 
and  direct  numerical  simulation  in  an  integrated  approach,  with  a  view  towards  the 
reduction  of  turbulent  sound  production  and  turbulent  drag. 

For  this  purpose,  in  a  first  step  that  is  targeted  at  identifying  interesting  domains 
in  the  space  of  parameters  describing  properties  of  a  compliant  wall  coating,  we  are 
developing  low-dimensional  models  based  on  Galerkin  projection  of  the  Navier-Stokes 
equations  onto  systems  of  eigenfunctions  obtained  via  Proper  Orthogonal  Decomposi¬ 
tion.  Because  of  the  relatively  small  effort  involved  in  simulating  and  analyzing  such 
models,  this  will  allow  us  to  scan  large  regions  of  parameter  space,  allowing  us  to  find 
regions  that  lead  to  a  reduction  of  turbulent  drag  and  turbulent  sound  production. 

Among  the  ultimate  goals  of  this  project  are  thus,  first,  to  obtain  a  fundamental 
understanding  of  flow-structure  interaction  phenomena  for  the  case  of  the  compliant- 
wall/turbulence  interaction,  and  second,  to  use  this  understanding  to  enhance  the  flight 
performance  of  air  vehicles  by  increasing  their  lift-to-drag  ratio. 


1  Statement  of  Objectives 

In  this  study,  research  is  conducted  into  the  issue  of  the  interaction  between  a  turbulent 
boundary  layer  and  an  adjacent  compliant  surface.  This  work  has  as  one  of  its  main  objec¬ 
tives  an  improved  understanding  of  the  dynamical  mechanisms  at  work  in  such  a  situation. 
Ultimately,  this  understanding  then  will  be  used  to  manipulate  the  turbulent  flow  such  that 
some  of  its  characteristics  are  altered  in  a  desirable  way.  In  this  respect,  we  are  mostly 
interested  in  two  goals.  We  want  to  be  able  to  reduce  the  turbulent  drag  generated  by  such 
boundary  layers,  and  we  also  want  to  be  able  to  attenuate  the  noise  that  is  generated  by  the 
turbulent  flow. 

In  order  to  achieve  this,  we  will  try  to  find  regions  in  the  space  of  parameters  describing 
mechanical  properties  of  compliant  coatings  (stiffness,  damping)  within  which  the  interaction 
between  the  wall  and  the  turbulent  flow  is  such  that  the  above  goals  can  be  achieved. 

The  research  proposed  here  consists  of  two  main  components.  To  be  able  to  determine 
regions  in  the  parameter  space  of  the  wall  coating  that  are  of  interest,  we  will  construct 
low-dimensional  models  based  on  Galerkin  projections  onto  Karhunen-Loeve  eigenfunctions. 
These  models  can  only  approximately  describe  the  dynamics  of  the  turbulent  flow,  but 
because  of  their  low-dimensionality,  they  allow  us  to  look  at  large  regions  of  the  parameter 
space.  Once  promising  combinations  of  parameters  are  determined,  we  will  then  use  accurate 
direct  numerical  simulations  to  assess  the  interaction  between  turbulence  and  the  compliant 
wall  in  more  detail. 

2  Status  of  Effort 

Our  work  on  developing  and  incorporating  models  for  the  interaction  between  the  turbulent 
flow  and  the  compliant  wall  is  proceeding  as  planned.  We  have  developed  a  refined  low¬ 
dimensional  model,  that  removes  many  of  the  limitations  of  the  previous  one,  results  of  which 
were  reported  in  our  previous  report.  In  parallel,  we  have  been  working  on  a  computer  code 
for  the  direct  numerical  simulation  of  turbulence/compliant-wall  interaction.  Results  from 
numerical  simulations  performed  with  that  code  are  also  presented  below. 


3  Accomplishments 

3.1  Introduction 

We  start  this  section  with  a  discussion  of  the  approximations  and  assumptions  that  go  in  a 
particular  low-dimensional  model  that  is  described  below.  We  will  first  present  an  order-of- 
magnitude  analysis  of  the  boundary  condition,  and  then  develop  an  implementation  of  the 
boundary  condition  in  the  pressure  term  under  general  conditions.  It  will  be  shown  that  for 
our  case,  as  a  reasonable  approximation  one  can  model  a  wall  that  moves  in  the  wall- normal 
coordinate  only,  which  allows  us  to  deal  with  non-zero  velocity  boundary  conditions  for  the 
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wall-normal  and  streamwise  coordinate  only.  We  will  also  derive  the  coupling  coefficient 
between  the  wall  motion  and  the  fluid  mechanics  that  appears  as  an  additional  term  in  the 
equation  for  the  wall-normal  velocity  component.  The  net  result  of  this  analysis  suggests 
that  for  the  purposes  of  our  low-dimensional  model,  neglect  of  the  nonlinear  terms  in  the 
boundary  conditions  is  usually  quite  justifiable. 

3.2  Construction  of  a  Low-Dimensional  Model  for  the  Wall-Flow 
Interaction 

In  this  section  we  will  introduce  the  dynamical  model  we  have  developed  for  describing  the 
behavior  of  near-wall  turbulence  over  a  compliant  boundary.  The  dynamical  equations  are 
obtained  from: 

•  Galerkin  projection  of  the  Navier-Stokes  equations  onto  a  special  set  of  basis  functions. 
In  order  to  capture  as  much  of  the  dynamics  of  the  flow  in  a  set  of  modes  that  is  as 
small  as  possible,  we  are  using  Karhunen-Loeve  eigenfunctions  (derived  from  a  rigid- 
wall  turbulent  flow),  which  are  complemented  by  Stokes  eigenfunctions  in  order  to  be 
able  to  account  for  the  compliant-wall  boundary  condition. 

•  The  eissumption  that  the  compliant  wall  can  be  modeled  as  a  simple  damped  mass¬ 
spring  system.  This  means  that  the  type  of  wall-coating  we  have  in  mind  can  be 
envisioned  as  some  sort  of  a  thin  rubber  layer,  and  it  will  become  clear  below  that  we 
are  also  assuming  the  deformations  of  the  wall  to  be  small.  These  assumptions  are 
compatible  with  the  results  reported  in  [3].  The  deformation  of  the  wall  is  driven  by 
pressure  fluctuations  (for  the  wall  normal  deformations),  and  by  the  fluctuating  wall 
shear  for  tangential  deformations  of  the  wall. 

3.2.1  Flow  Equations 

For  these  preliminary  calculations,  we  are  constructing  a  model  that  is  analogous  to  the 
one  first  presented  by  Aubry  et  al.  [1],  which  is  also  described  in  much  greater  detail  in 
the  monograph  by  Holmes,  Lumley  &  Berkooz  [5].  This  means  in  particular  that  we  are 
dealing  with  streamwise-invariant  modes  only,  arguing  that  these  modes  should  dominate 
the  behavior  of  our  flow  just  as  they  have  been  demonstrated  to  do  in  the  rigid-wall  case. 
Thus,  we  are  going  to  represent  our  flow  fields  using  a  decomposition  of  the  form 

u(y,z,t)  =  I7(y)  +  5^a(^'’'=V^'’''n2/)e‘''"  +  c.c.,  (1) 

j,k 

where  y  and  2  are  the  spanwise  and  wall-normal  coordinates,  respectively,  k  is  the  spanwise 
wave  number,  and  j  denotes  the  order  of  the  basis  function  corresponding  to  a  given  wave 
number  k.  The  function  jg  jth  basis  function  corresponding  to  wave  number  k. 

In  the  model  described  below,  the  index  j  will  assume  the  values  of  1  or  2  only.  This 
means  that,  as  in  [1],  we  will  be  using  just  one  POD  mode  per  wave  number,  described 
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by  the  function  For  the  second  and  third  mode,  we  will  use  Stokes 

eigenfunctions  (see  below,  §3.2.5).  The  main  role  of  these  Stokes  functions  will  be  to  allow 
us  to  match  non-zero  velocity  boundary  conditions  at  the  wall,  which  are  necessary  for  a 
moving  w^all  (remember  that  the  are  POD  modes  from  a  rigid- wall  case,  w’hich  have 

=  0)  =  0).  The  Stokes  functions  are  chosen  such  that  they  are  orthogonal  to  the 
Karhunen-Loeve  modes,  so  that 

With  these  preparations,  the  equations  for  the  fixed  wall  (for  each  spanwise  w'avenum- 
ber  k.  and  for  each  “quantum  number’'  j)  are  of  the  form  [5] 


a 


(kk) 


+  (1  + 


+  ^2  ^(k\k-k'), 
k',p,q 


Q^(p.k')^{q,k-k>) 


n,k',p,q 


(2) 


Since  we  are  actually  interested  in  dealing  not  with  a  fixed  but  wdth  a  deformable  w^all,  the 
above  equations,  derived  for  a  fixed  wall,  are  not  sufficient  to  describe  that  situation,  and 
need  to  be  extended  by 

•  adding  equations  that  describe  the  new  boundary  conditions  for  a  moving  wall,  and 

•  taking  into  account  the  influence  of  the  pressure  term  at  the  moving  w^all. 


3.2.2  Boundary  Conditions 

Our  coordinate  system  is  such  that  the  undisturbed  compliant  surface  occupies  the  plane  y  = 
0.  Let  ^{z)  =  [^1.  ^2- be  the  displacement  vector  of  the  compliant  surface  (remember 
that  w'e  assume  streamwise  invariance).  The  correct  set  of  boundary  conditions  is  given  by 

Ui{x  +  ^,t)  =  ii,  (3) 

where  u  is  the  instantaneous  fluid  velocity,  and  ^  denotes  the  time  derivative  of  the  w^all 
displacement  vector.  For  our  model,  we  will  aissume  that  the  surface  displacements  are  small, 
more  specifically,  we  assume  that  |£''’|  =  ^  5,  so  that  the  velocity  profiles  are  linearL 

Under  this  assumption,  the  boundary  condition  can  be  linearized,  giving 

Ui{{i,z,t) ^  ii{z,t).  (4) 

If  w-e  express  the  instantaneous  fluid  velocity  in  terms  of  the  fluctuating  components  Ui{x,t) 
and  the  mean  U  =  [f/i(i/),  0, 0], 

Ui  ~  Ui  -f-  tXj,  (5) 

'Note  that  this  magnitude  of  ^  coincides  with  the  findings  in  [3]. 
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we  can  write  (4)  as 

«,(0,  2,  t)  +  [/^J(0)e,•(2,  t)  +  Ui,,(0,  2,  t)  =  U-,  t).  (6) 

where  we  have  used  Ui{Q)  =  0.  Let  us  now  assume  that  Ui(0,  z,  t)  and  ^i(z.  t)  are  both  small 
quantities  in  a  sense  yet  to  be  defined.  Let  us  consider  the  velocity  gradient  tensor  at  y  =  0. 
u,j(0.  z.t).  We  may  write 

/  0  Wi,2  «1.3  \ 

0  ^  (7) 

^  0  tt3,2  ! 

where  we  have  indicated  by  an  underbrace  terms  that  will  vanish  for  vanishing  wall  defoi- 
mation  ^  >  0.  Note  that  for  the  rigid-wall  case  ^  =  0,  we  have  ^2,2  =  0  from  continuity, 

and  Ui,3  =  0  because  the  velocity  vanishes  identically  at  the  wall.  We  now  proceed  to  elim¬ 
inate  all  those  terms  in  (4)  that  are  quadratic  in  quantities  that  vanish  when  ^  =  0.  The 
precise  magnitude  of  the  errors  that  this  step  incurs  will  be  discussed  below.  The  boundary 
conditions  with  the  quadratic  terms  eliminated  will  be  referred  to  as  the  "ultra-linearized 
boundary  conditions, 

u\  +  f/1,2^2  +  t^l,2^2  =  ?1) 

U2  “  ^2,  (.8) 

U3  ■+-  'U3,2^2  =  Cs- 

Note  that  the  above  set  of  boundary  conditions  still  contains  two  terms,  Ui,2^2  and  U3.2C21 
that  are  nonlinear  in  the  dependent  variables  of  our  problem.  Because  of  the  complications 
associated  with  this  nonlinearity,  we  want  to  drop  these  terms,  too.  We  can  estimate  their 
magnitude  from  experimental  data  in  Townsend  [7].  There  we  find  that  at  y'^  =  5,  where 
y+  zz:  UrV/i^.  Ui/Ur  -  1.54,  while  U^/Ur  =  0.306.  The  behavior  of  both  is  linear  between 
y'*'  =  5  and  the  wall.  We  may  thus  estimate 

■Ui  2  ~  0.308 — 

’  u 

U3  2  ~  0.0612-^  (9) 

u 


for  values  of  y"*”  between  0  and  5.  Since  we  have  U\^2  ~  see  that  these  derivatives 

(at  least,  uj  2)  are  not  substantially  smaller  than  the  mean  velocity  derivative.  However,  if 
we  accept  an  accuracy  of  30%,  we  can  write 


til  +  ^^1,2^2 


=  6 
=  iz- 


We  can  now  estimate  the  orders  of  magnitude  of  the  three  components. 


-6,  «3 
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and  we  can  also  accurately  estimate  the  magnitude  of  the  omitted  components  (the  ones  that 
were  quadratic  in  terms  which  vanish  when  ^  vanishes).  A  single  example  will  be  sufRcient: 

^  »  &*:3,  (12) 

Ui 

and  the  others  are  similar,  ks  is  the  wavenumber  in  the  2-direction.  Thus,  these  terms  will 
be  small  if  the  displacement  of  the  surface  is  small  relative  to  the  wavelength  of  the  surface 
disturbance.  This  will  almost  always  be  the  case.  Finally,  let  us  consider  the  terms  that  do 
not  vanish  when  ^  vanishes.  We  obtain 

0.308— ,  «  0.0612|^.  (13) 

Ui  4i  “3  « 

We  have  thus  shown  that  by  neglecting  relative  to  1/1,2^21  we  make  a  maximum 

error  of  about  30%.  Henceforth,  for  the  low-dimensional  model,  we  will  assume  (10).  because 
of  its  linearity.  The  relatively  large  potential  error  that  is  involved  in  this  assumption  is 
one  of  the  reasons  why  this  project  will  also  include  direct  numerical  simulations  of  the 
turbulent  .flow  over  a  compliant  boundary,  in  order  to  be  able  to  assess  the  effect  of  the 
various  simplifications  involved  in  the  modeling  process.  Of  course,  for  the  direct  numerical 
simulations  that  we  are  planning,  the  above  two  terms  will  be  kept.  Since  we  are  using  a 
dynamical  coordinate  transformation  technique  for  this  simulation  (see  below),  there  is  no 
difficulty  in  treating  the  boundary  conditions  exactly,  so  the  above  approximations  are  not 
necessary  in  that  case. 


3.2.3  The  Pressure  Term 

In  the  dynamical  systems  model  for  the  flow  as  given  by  (2)  above,  the  pressure  term  in 
the  Navier-Stokes  equations  has  been  neglected.  Formally,  after  Galerkin  projection,  the 
pressure  term  results  in  an  expression  of  the  form 

grad  (p)^  =  J  •  grad  (p)  dec,  (14) 

where  is  one  of  the  three-dimensional,  vector-valued  basis  functions  being  used  in  the 

projection  (in  our  case,  z)  —  <p>b’^l(p)e’*^^),  and  the  integral  is  over  the  flow  domain 

of  interest.  Using  the  integral  theorems  and  the  fact  that  the  2)  are  solenoidal,  this 

can  be  written  as  a  surface  integral 


where  n  is  the  outward  normal  on  the  boundary  S  of  the  domain.  Using  a  Fourier  repre¬ 
sentation  of  the  pressure,  p{y,z)  =  '^i^p{k,y)e^'^^ ,  we  find  that  the  only  remaining  terms 


,grad(p)')  =  -pik,0)ip^/*{0)  +  p{k,ym)^^2''*iyn 
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where  ym  is  the  height  of  our  domain. 

In  [1.  5]  it  is  shown  that  the  term  from  the  upper  boundary  mainly  acts  as  a  random 
disturbance  to  keep  the  solution  at  a  finite  average  distance  in  phase  space  from  the  attracting 
heteroclinic  orbit,  which  is  responsible  for  the  bursting  events  in  the  near-wall  turbulent  flow. 
It  is  demonstrated  there  that  this  pressure  term,  however,  does  not  fundamentally  alter  the 
dynamics  of  the  system.  In  addition,  since  we  are  interested  only  in  comparing  the  behavior 
of  the  compliant-wall  system  to  the  one  with  the  fixed  wall,  omission  of  the  last  term  in  (16) 
is  justified.  However,  the  first  term  on  the  right-hand  side  of  (16)  is  zero  only  for  the  fixed 
wall,  but  will  in  general  be  non-zero  for  a  moving  wall.  We  therefore  want  to  add  this  term 
to  our  dynamical  equations,  so  that  we  end  up  with 


d 


O'.fc) 


Y.  ('■i-r""' + (1 + iTTow-r) 


n,k',p,q  ^ 


(17) 


3.2.4  Equations  for  the  Dynamics  of  the  Wall 

We  model  the  compliant  surface  in  the  Fourier  space  spanned  by  the  spanwise  wavenumbers  k 
as  a  linear  damped  mass-spring  system  driven  by  the  fluid  stress  at  the  wall.  The  stiffness 
and  damping  parameters  may  be  different  for  different  spanwise  wavenumbers  k.  In  the 
streamwise  and  spanwise  directions,  the  deformations  of  the  wall  are  driven  by  the  shear 
stress  which  is  easily  accessible  in  the  models  as  a  linear  term  in  the  model  coefficients. 


Mflf  -f  +  kU'I  =  r^Au,aii  = 
+  Dlh  +  =  r^Uii  = 


n 

n 


(18) 

(19) 


The  wall-normal  deformations  of  the  compliant  surface  are  driven  by  the  pressure  fluctuations 
at  the  boundary. 

=  p(0,  k)  (20) 

We  employ  capital  letters  to  describe  the  properties  of  the  compliant  surface  for  clarity; 
the  lower  indices  allow  for  the  elastic  properties  of  the  wall  to  be  different  for  different 
coordinates. 


3.2.5  Basis  Functions 

We  choose  to  base  our  low-dimensional  model  for  the  turbulent  boundary  layer  over  a  compli¬ 
ant  surface  on  the  eigenfunctions  of  the  rigid-wall  boundary  layer.  Since  these  eigenfunctions 
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have  =  0  at  the  wall  and  thus  do  not  allow  us  to  satisfy  the  boundary  conditions  (10). 
we  need  to  introduce  additional  eigenfunctions  to  account  for  the  motion  of  the  compliant 
surface.  In  the  absence  of  surface  compliance,  our  models  will  revert  to  the  rigid-wall  case. 
Lacking  an  experimental  or  computational  database  on  the  turbulent  flow  ov^er  a  compliant 
surface,  we  are  reverting  to  derive  our  additional  eigenfunctions  in  an  ad  hoc  fashion.  We 
choose  additional  eigenfunctions  as  solutions  of  the  Stokes  equation  with  periodic  motion  of 
the  wall. 

W'e  will  first  approach  the  Stokes  problem  with  boundary  conditions  corresponding  to 
the  oscillation  of  the  wall  in  the  streamwise  direction.  Since  our  low-dimensional  model  for 
the  rigid-wall  boundary  layer  is  built  upon  a  Fourier  decomposition  of  the  velocity  field,  we 
will  use  the  Stokes  equation  to  generate  an  eigenfunction  for  each  wavenumber  k\ 


A-ik|»  =  0, 

u  ot 


,2  1  ^ 

at 


u 


{k,y,t)  =  0, 


(21) 


2  ^ 

and  we  will  then  proceed  to  assign  these  solutions  to  our  basis  functions  ^Pi  .  The  velocity 
field  is  subject  to  the  boundary  conditions 


u{k,0,z)  =  cos{Pt), 

u{k,  y  ^  oo,t)  — >  0. 


(22) 

(23) 


In  order  to  satisfy  the  oscillating  boundary  conditions,  the  velocity  field  u{y,t)  can  be 
represented  as  a  superposition  of  left-  and  right-running  traveling  waves. 

u{k,y,t)  =  u^{k,y)e'i^^  +  u-{k,y)e-^^k  (24) 

We  may  sol\-e  separately  for  u~{k,  y)  and  u'''{k,  y)  from  the  governing  equation  and  boundary 
conditions. 


u+{k,y)  =  +  (25) 

u-{k,y)  =  +  (26) 


In  order  for  u(y,  t)  to  decay  at  infinity,  the  first  term  in  both  u'^  and  u~  must  disappear, 
so  that  A'^  =  A~  =  0.  The  boundary  condition  at  the  wall  may  be  satisfied  by  and  u~ 
together. 


u{k,y  =  0,t) 


u+{k,  y  =  0,  t)e‘^‘  +  u-{k,  y  =  0,  t)e--^‘ 
cos{Pt) 
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(27) 

(28) 
(29) 


The  choice  \  satisfies  the  boundary  conditions  and  results  in  the  following  form 

for  u{y.  t): 

u{k,y,t)  =  (30) 

A  ^ 

(31) 


We  choose  the  eigenfunction  which  represents  streamwise  compliance  from  the  initial  value 
of  u{k.  yA).  i.e.  ■^f^\y)  =  u{k,yA  =  0),  <pP^(2/)  =  =  0.  This  choice  for  the  form 

of  the  eigenfunction  will  permit  the  boundary  conditions  on  the  streamwise  component  of 
velocity  to  be  satisfied. 

We  would  ideally  prefer  each  additional  eigenfunction  to  have  only  a  single  component 
with  a  non-zero  boundary  condition.  This  will  not  always  be  possible  when  we  account  for 
the  motion  of  the  wall  in  all  three  directions,  and  the  eigenfunctions  are  orthogonalized. 
But  in  the  current  case  of  strictly  wall-normal  motions  and  using  the  "ultra-linearized 
boundary  conditions  (10),  the  additional  eigenfunctions  may  , be  formulated  in  this  fashion 
with  =  0)  ^  0  and  =  0)  7^  0.  All  other  components  are  zero  at  the  wall. 

In  fact,  the  wall-normal  and  spanwise  components  of  and  the  streamwise  component  of 
are  strictly  zero. 

We  now  determine  the  additional  eigenfunctions  for  the  wall-normal  motion  of  the  wall 
in  a  similar  fashion,  using  the  streamfunction  formulation  of  the  Stokes  equation: 


aj/2 


^'(A:,y,t)  =  0 


(32) 


As  before,  we  superpose  left-  and  right-running  wave  solutions  to  this  equation  '^{k.yA)  = 
e'‘^'^'*'(/c.  y)  -^ e"''^^ y),  and  we  apply  the  boundary  conditions  at  infinity  to  and  'I'" 
to  eliminate  the  terms  which  grow  exponentially.  Next,  we  enforce  the  boundary  condition  011 
the  spanwise  velocity  d'^{y  =  0;k)/dy  =  0.  The  two  wave  solutions  are  combined  to  satisfy 
the  time-varying  boundary  conditions  at  the  wall,  ik'^{y  =  0,/c)  =  cos(/?t).  This  procedure 
allows  the  formation  of  a  divergence-free  additional  eigenfunction  which  in  turn  allows  the 
low-dimensional  models  to  satisfy  the  boundary  conditions  at  the  compliant  surface.  The 
resulting  basis  functions  are  shown  in  Fig.  1. 

After  the  additional  eigenfunctions  are  obtained  as  solutions  of  Stokes  problems,  the 
complete  set  of  basis  functions  is  orthogonalized  using  a  Gram-Schmidt  procedure.  This 
gives 


(2,fc) 


(  (2,k) 

(  (l,fc)  ’ 


0, 

0, 


(33) 
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and 


(3,fc1 

'A 


V^2 


'r'i 


0. 


,J3,fc) 

V2 


(3,fc) 


(  (3,fe)  (l,fc)\  ,  (  (3,fc)  ,  Jl.fc)^ 

\<f2  )  +  [f3  .^3  j 

'3  /  (S.fc)  (3,fc)\  ^  /,,(3,fc)  ,J3M)\  ■ 

(^V^2  J  +  \^V^3  )V^3  j 


(34) 


The  new  basis  functions  are  finally  normalized  so  that  =  1. 


3.2.6  Six-Mode  Model  with  Wall-Normal  Compliance 

The  relations  (10),  (17),  and  (18-20)  constitute  a  complete  set  of  evolution  equations  for  the 
coefficients  of  the  eigenfunctions  as  well  as  evolution  equations  for  the  motion  of  the  wall. 
In  the  following,  we  will  further  simplify  this  model  by  assuming  that  the  only  significant 
motion  of  the  wall  is  in  the  wall-normal  direction,  which  means  that  we  will  drop  (18) 
and  (19).  In  this  case  we  need  only  two  additional  basis  functions  (the  Stokes  modes  given 
above)  to  satisfy  the  wall-normal  and  streamwise  boundary  conditions.  This  way  we  end  up 
with  a  model  that  includes  one  POD  mode  from  the  fixed-wall  case,  one  Stokes  mode  for 
the  wall-normal  and  streamwise  coordinate  boundary  conditions,  respectively,  and  one  mode 
for  the  motion  of  the  wall,  which  is  later  eliminated  by  using  the  coupling  between  wall  and 
fluid  motion.  Each  of  the  three  wall-normal  modes  is  then  multiplied  by  one  of  two  Fourier 
modes,  which  brings  the  total  to  six  complex  modes  that  make  up  this  model. 

We  choose  to  determine  the  coefficients  for  the  additional  eigenfunctions  fiom  the  bound¬ 
ary  conditions  at  the  wall.  The  pressure  may  then  be  determined  from  the  evolution  equation 
for  the  additional  eigenfunction  representing  the  wall-normal  motion  of  the  wall 


p{k.y  =  0) 


=  0) 


t) 


c 


(3) 

(k',k-k') 


PQ 


k',p,q 


(35) 


The  cubic  terms  disappear  because  this  additional  eigenfunction  is  defined  to  have  no  stream- 
wise  component,  and  there  is  no  streamwise  variation  in  our  model.  When  this  expression  is 
substituted  into  the  equation  for  the  wall-normal  motion  of  the  wall  (20),  the  time  deriva¬ 
tive  term  merges  into  the  resulting  in  a  new  effective  mass  of  the  wall  M^-  In  fact, 

^2  =  since  is  the  only  eigenfunction  with  a  non-zero  vertical  velocity  at  the 
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wall.  Normalizing  with  this  new  mass  using 


M  = 

M+p/(4'>2r)iu., 

(36) 

P  = 

p/(M4f  U), 

(37) 

D  = 

D/M, 

(38) 

k  = 

K/M, 

(39) 

results  in  the  following  low-dimensional  model  for  the  flow  over  a  compliant  surface; 


^(p,fc')a(9.fc-fc') 


r.k^pyQ 


I2  +  + -^2  4  ~ 

.  p 

(41) 

k' ,p,q 

ck  dU\ 

S2  dy  l!/=0 

(42) 

^(3, A:) 

(43) 

3.3  Results 

In  the  following,  we  will  report  the  findings  we  obtained  from  an  investigation  of  the  above 
six-mode  model,  where  we  have  used  two  spanwise  Fourier  modes  k  for  each  of  the  basis 
functions.  As  a  result,  we  are  looking  at  the  behavior  of  a  model  described  by  a  system  of 
ODEs  of  order  12. 

We  started  our  investigation  of  this  model  by  systematically  varying  p,  D  and  K  as  our 
input  parameters,  and  monitoring  the  resulting  bursting  frequency  as  our  output  parameter. 
The  role  of  p  is  to  determine  the  relative  strength  of  the  wall  forcing.  Increasing  p  increases 
the  amplitude  of  the  wall  motion,  and  vice  versa.  We  tried  values  of  p  between  p  —  0.01 
and  p  =  1,  within  which  range  our  findings  do  not  change  qualitatively.  For  values  of  p 
significantly  smaller  than  that,  the  wall  motion  is  too  small  to  have  a  noticeable  effect,  and 
for  values  of  p  larger  than  that  range,  the  dynamics  of  our  system  changes  qualitatively,  so 
that  the  appropriateness  of  the  model  becomes  questionable. 


We  found  that  the  damping  constant  D  has  relatively  little  effect  (the  reason  for  this 
will  become  clear  below),  whereas  the  stiffness  K  turned  out  to  be  the  most  important 
parameter.  Varying  this  stiffness  parameter  we  found  that  the  bursting  frequency  could  be 
reduced  by  as  much  as  30%.  We  note  that  one  would  expect  this  to  correspond  to  a  drag 
reduction  of  the  same  order,  although  there  is  a  possibility  that  the  intensity  of  the  bursts 
changes  through  the  influence  of  the  compliant  wall,  which  could  modify  the  actual  value  of 
drag  reduction  that  is  realized.  While  we  do  not  see  a  change  in  the  quality  of  the  bursts 
in  our  model,  this  might  be  an  artifact  of  the  modeling  process,  which  is  another  reason  for 
us  to  check  the  findings  from  our  model  through  a  direct  simulation  that  we  are  currently 
working  on. 

Fig.  2  shows  a  comparison  of  the  phase-space  portraits  for  the  cases  of  a  rigid  and  a 
compliant  wall,  respectively.  The  parameters  we  have  chosen  for  the  compliant  wall  are  the 
ones  that  give  minimum  bursting  frequency.  The  comparison  of  the  two  cases  clearly  shows 
that  the  qualitative  dynamics  of  the  flow  is  hardly  changed  at  all  by  the  presence  of  the 
compliant  wall.  This  is  important  because  it  gives  credence  to  our  assumption  that  we  can 
still  describe  the  velocity  distribution  of  the  compliant-wall  flow  by  using  the  most  energetic 
eigenmode  from  a  fixed-wall  case. 


Figure  2;  Phase-space  projection  onto  the  plane  spanned  by  the  real  parts  of  the  coefficients 
of  the  POD  modes  for  the  first  and  second  Fourier  mode,  respectively,  a),  fixed  wall.  b). 
compliant  wall.  Parameters  where  d  =  0.1,  p  =  0.1,  k  =  500. 


Fig.  3  shows  a  comparison  of  the  behavior  in  time  of  the  coefficient  of  the  most  energetic 
POD  mode.  These  time  traces  exemplify  the  intermittent  nature  of  the  flow  dynamics:  the 
flow  stays  calm  for  extended  periods  of  time  (the  plateaus  in  Fig.  3),  corresponding  to  the 
rolls  being  stationary,  or  very  slowly  drifting  from  side  to  side.  Then,  during  a  burst  (the 
steep  gradients  in  Fig.  3),  the  rolls  begin  to  move  violently,  break  down,  and  are  recreated 
with  a  phase  shift  of  a  half  wavelength. 
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Most  significantly,  however,  the  figure  shows  that  for  the  compliant-wall  case,  the  quies¬ 
cent  periods  are  longer,  which  corresponds  to  a  reduced  bursting  frequency.  Quantitatively, 
it  turns  out  that  in  this  case  the  bursting  frequency  has  been  reduced  by  about  307c.  We 
should  alert  the  reader  to  the  fact  that  the  increase  in  inter-burst  time  for  the  compliant-wall 
case  is  a  transient  effect  for  the  numerical  simulations  shown  in  the  figures.  As  discussed 
in  '1]  (see  also  [5]  for  a  more  detailed  discussion),  for  a  model  without  external  disturbances, 
the  inter-burst  times  should  become  longer  and  longer  as  the  system  gets  closer  to  the  hete¬ 
roclinic  cycle.  However,  due  to  finite  truncation  errors,  any  numerical  simulation  will  show  a 
convergence  of  the  inter-burst  times  towards  a  finite  value.  The  ultimate  value  that  one  finds 
in  a  numerical  solution  is  a  function  of  the  numerical  accuracy  of  the  integration  scheme  that 
was  used  and  has  no  physical  meaning.  In  a  real  flow,  on  the  other  hand,  there  will  be  pres¬ 
sure  disturbances  acting  on  the  upper  boundary  of  the  near-wall  region  we  are  considering 
here.  In  that  case,  inter-burst  times  depend  on  the  magnitude  of  these  pressure  disturbances 
as  well  as  on  how  strongly  attracting  the  heteroclinic  cycle  is.  If  we  have  a  heteroclinic 
cycle  that  is  more  strongly  attracting,  then,  on  average,  the  solution  will  be  closer  to  the 
heteroclinic  cycle,  which  results  in  longer  inter-burst  times.  This  is  precisely  the  scenario 
that  is  indicated  by  the  longer  inter-burst  times  we  observed  during  the  transients  in  our 
numerical  experiments.  The  prolonged  inter-burst  periods  in  the  transients  demonstrate 
that  the  heteroclinic  cycle  for  the  compliant  wall  model  is  more  attracting,  and  hence,  if 
numerical  accuracy  were  not  a  factor,  inter-burst  times  would  consistently  be  higher  for  the 
compliant-wall  case. 

Finally,  we  should  note  that  within  our  numerical  accuracy  the  location  of  the  heteroclinic 
cycles  is  identical  for  the  rigid-wall  and  the  compliant-wall  case.  For  Fig.  3  we  have  started 
from  identical  initial  conditions  for  both  the  fixed-wall  and  the  compliant-wall  case,  which 
means  that  the  results  shown  in  that  figure  do  indeed  demonstrate  that  the  heteroclinic  cycle 
for  the  compliant-wall  case  i.s  more  strongly  attracting  than  the  one  for  the  fixed  wall.  In 
further  numerical  experiments  with  different  initial  conditions  (not  shown  here)  we  found 
qualitatively  the  same  behavior. 

Our  last  figure  shows  a  representation  of  the  temporal  dynamics  of  the  wall  motion  itself. 
W'e  can  sec  that  the  frequency  of  the  wall  motion  is  about  a  factor  of  50  higher  than  the 
bursting  frequency.  In  addition,  the  reader  should  take  note  of  the  extremely  small  amplitude 
of  the  wall  fluctuations  as  shown  in  Fig.  4.  This  observation  validates  our  assumption  of 
small  wall  deformation  we  made  in  deriving  our  model,  and  it  also  demonstrates  why  the 
value  of  the  damping  constant  in  our  model  has  little  influence  on  the  dynamics  of  the  flow. 
At  these  tiny  fluctuation  levels,  the  energy  that  is  dissipated  in  the  compliant  wall  is  many 
orders  of  magnitude  lower  that  the  energy  dissipation  of  the  flow,  so  one  would  expect  to 
see  next  to  no  influence  of  energy  dissipation  in  the  wall  on  the  flow  dynamics. 

3.3.1  Direct  Numerical  Simulation 

The  problem  we  are  considering  is  basically  a  fluid-solid  problem.  On  the  fluid  side,  we 
have  the  incompressible  NS  equation  as  the  model  for  fluid  motion.  On  the  solid  side,  we 
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Time 

Figure  3:  Comparison  of  the  time  behavior  of  the  real  part  of  the  coefficient  of  the  POD  mode 
for  the  second  Fourier  mode.  Solid,  compliant  wall,  dashed,  fixed  wall.  Same  parameters  for 
the  compliant  wall  as  above. 


use  a  spring-supported  elastic  plate  to  model  the  compliant  wall.  The  equation  governing 
the  w'all-normal  motion  nondimensionalized  by  the  average  half  channel  height  h  and  the 
laminar  centerline  velocity  uq  reads 


dri 


(d>r,  ,  d'n  a‘r,  \ 
dxj  dx\xl) 


+  C3TI  —  Cx 


d'^T) 


n  -  _ 
^^dxl 


(44) 


where  r]  is  the  wall-normal  displacement,  is  the  pressure  disturbance  on  the  wall  and  Co- 
C].  C2.  C3,  and  Cx{Cz)  are  the  nondimensional  wall  properties  defined  as 


Co  = 


Pm^ 


c,= 


1 


ph 

Ksh^ 


r  -  —— 

^  Repu’ 


Re^  pu^ 


Cx  = 


1  Txh 


0.= 


J _ ^ 

Re^  phu'^ '' 

1  Tzh 
Re^  pv"^  ’ 


(45) 


Re^  pu^ 

The  wall  parameters  are:  the  plate  density  pm  and  thickness  b,  the  wall  damping  coefficient 
d.  the  flexural  rigidity  of  the  plate  B,  the  streamwise  (spanwise)  tension  per  unit  width 
Tx{Tz)  and  the  spring  stiffness  Ke]  P  is  the  density  of  the  fluid  and  u  is  its  kinematic 
viscosity.  The  turbulent  channel  flow  we  simulate  is  assumed  periodic  in  the  streamwise(Ti 
x)  and  spanwise(x3  or  z)  directions  and  has  a  rigid  upper  and  a  compliant  lower  wall  that 


or 


are  separated  by  an  average  distance  of  2h  in  the  wall-normal(a:2  or  y)  direction.  Re  is  its 


Reynolds  number  based  on  h  and  uq. 

The  main  difficulty  of  the  simulation  is  to  handle  the  boundary  conditions  at  the  moving 
compliant  wall.  use  a  time-varying  coordinate  transformation  to  eliminate  the  deforma¬ 
tion  of  the  compliant  boundary  in  the  computational  domain.  The  continuity  equation  in 
the  computational  coordinate  system  becomes 


u 


i 


=  0, 


(46) 


and  the  momentum  equation  in  the  rotational  form  becomes 

dn'  dT-^  du^  ,  „  dp  ^  1 

ar  =  -— »  to- + 


(47) 


where  u’  is  a  contravariant  velocity  component,  the  semicolon  denotes  covariant  differentia¬ 
tion,  T~^  is  the  inverse  coordinate  transformation  function  in  the  wall-normal  direction,  and 
g’^  is  the  contravariant  metric  tensor  of  the  transformation.  The  equation  for  pressure  p  is 
obtained  by  taking  a  covariant  differentiation  with  respect  to  Xi  on  the  momentum  equation 
and  can  be  written  as 


dT-^  du' 
dt  dx2 


4“  Vp g 


(48) 


The  boundary  conditions  for  velocity  must  also  be  transformed  into  the  new  coordi¬ 
nate  system.  The  boundary  conditions  for  pressure  can  then  be  determined  indirectly  from 
the  incompressibility  condition  via  an  influence  matrix  technique.  A  Fourier-Galerkin  and 
Chebyshev-Tau  pseudospectral  method  for  spatial  discretization  and  a  three-sub-step  RK 
method  for  time  advancement  are  used  for  solving  the  equations. 
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Results  The  computational  domain  sizes  in  the  streamwise  and  spanwise  directions  are 
chosen  to  be  tt  and  respectively.  The  computation  is  carried  out  using  32  x  65  x  32  (in 
X,  y,  z)  grid  points  for  a  Reynolds  number  of  3300.  In  addition,  as  a  reference  case  for 
comparisons  we  have  simulated  the  turbulent  flow  in  a  canonical  channel  composed  of  two 
rigid  walls  with  the  same  flow  parameters  as  those  for  the  compliant  wall.  For  both  cases 
the  mean  pressure  gradient  is  calculated  at  each  time  step  to  enforce  constant  mass  flux  in 
the  ctiannel. 

An  important  issue  in  the  compliant-wall  case  is  how  to  choose  the  compliant  wall  param¬ 
eters.  We  currently  choose  them  following  [Semenov  (1991)],  while  satisfying  the  numerical 
requirements  of  stability,  the  accuracy  of  the  representation  of  the  physical  fluid-solid  inter¬ 
action.  and  an  acceptable  numerical  time  step.  The  damping  factor  of  the  compliant  wall 
plays  a  key  role  in  the  simulations.  It  should  be  large  enough  to  inhibit  a  dynamic  instability 
of  the  interaction  between  the  compliant  wall  and  the  turbulent  flow.  On  the  other  hand, 
the  compliant  wall  should  move  intensively  enough  so  that  it  physically  interacts  with  the 
turbulence  rather  than  just  generating  some  numerical  noise  in  the  computational  system. 
These  requirements  lead  to  values  for  the  wall  parameters  that  are  not  very  close  to  those 
of  materials  that  others  have  used  in  their  experiments,  because  the  small  pressure  distur¬ 
bance  due  to  the  low  Reynolds  number  of  our  DNS  cannot  excite  a  sufficient  motion  of  an 
experimental  compliant  wall.  Nevertheless  we  believe  to  observe  similar  physics.  The  third 
numerical  requirement  is  that  the  natural  frequencies  corresponding  to  pressure  disturbances 
of  high  wave  numbers  have  to  be  low  enough  to  guarantee  an  acceptable  time  step  in  the 
numerical  simulation.  In  a  monoharmonic  approximation  for  the  compliant  wall  motion,  the 
forcing  pressure  disturbance  is  given  in  the  harmonic  form 

=  (49) 

The  nvaural  circular  frequency  of  the  free  vibration  /  and  the  time  constant  of  damping  C 
for  the  compliant  wall  are  then  found  to  be 


0  (no  free  vibration) 

4  1/?^  - 


if  /?  >  2a;, 
if  /3  <  2a;, 


(50) 


0-{0" 


where  (3  and  lo  are  defined  as 


if  /?  >  2a;, 
if  /?  <  2a;, 


(51) 


=  2-  +  kt  +  2klki)  +  c,kl  +  c,ki  +  c,] . 


(52) 


The  steady  amplitude  (in  wall  units)  of  the  compliant  wall  displacement  r'*'  and  the  normal 
velocity  of  the  wall  FJ"  are 


^  _  KdVwR^r 
CoLO 


(53) 
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Figure  5;  The  results  of  the  monoharmonic  analysis  for  the  compliant  wall  motion  for  Re  — 
3300.  Rcr  =  140,  =  0.01,  with  k2  =  kl  +  kl. 


y+  _  ^dPw(^R^ 

CouRer  ■ 

Here.  Re-r  is  the  turbulent  Reynolds  number  on  the  compliant  wall,  and  K^.  is  the  dynamic 
coefficient  which  can  be  calculated  from 

(55) 

Figure  5  shows  the  results  of  a  monoharmonic  analysis  for  the  compliant  wall  motion.  These 
and  the  later  simulation  results  indicate  that  the  compliant  wall  we  used  in  the  simulation 
satisfies  the  numerical  requirements.  The  flow  and  the  compliant  wall  parameters  in  the 
simulation  are  given  in  Table  1. 

Figure  6  shows  the  history  of  the  total  drag  (represented  by  the  mean  streamwise  pres¬ 
sure  gradient)  and  the  skin  friction  coefficients  for  the  two  cases.  The  skin  friction  on  the 
compliant  wall  is  12%  less  than  that  on  a  rigid  wall  in  the  reference  case,  while  they  are 
about  the  same  on  all  rigid  walls  in  both  cases. The  turbulent  Reynolds  number  is  138  based 


Kd  = 


1  - 


cr‘ 


-H 
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p  1000  kg 

h  0.005  m 

u  lx  10“® 

Re  3300 


800  kg  m“^ 
h  0.003  m 

d  2000  kg 

B  3  X  10“®  kg  m"^  s“^ 

Tx  0.001  kg  s-2 

T,  0.0001  kg  s-2 

Ke  10000  kg  m"^  s“^ 


Table  1:  Values  of  the  flow  and  the  compliant  wall  parameters. 
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Figure  6;  Time  traces  of  (a)  total  drag,  (b)  skin  friction  on  the  compliant  wall  in  the 
compliant  wall  case  and  a  rigid  wall  in  the  reference  case  and  (c)  skin  friction  on  the  rigid 
wall  in  the  compliant  wall  case  and  a  rigid  wall  in  the  reference  case,  solid  line:  the  reference 
case,  dashed  line:  the  compliant  wall  case. 
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Figure  7:  Comparisons  of  (a)  turbulent  kinetic  energy,  (b)  turbulent  kinetic  energy  produc¬ 
tion,  (c)  RMS  streamwise  vorticity  and  (d)  Reynolds  shear  stress.  The  average  position  of 
the  compliant  wall  is  at  y  =  —  1.  solid  line:  the  reference  case,  dashed  line:  the  compliant 
wall  case 


on  the  compliant  wall  friction  velocity  and  147  on  a  rigid  wall.  Consistent  with  these  obser¬ 
vations.  we  achieve  a  total  drag  reduction  of  6%  in  a  channel  with  one  compliant  wall.  The 
amplitude  of  the  compliant  wall  motion  is  so  small  that  most  of  the  present  drag  reduction 
is  due  to  the  wall-normal  velocity  of  the  compliant  wall  rather  than  the  wall  displacement 
and  the  total  drag  is  almost  completely  from  skin  friction.  It  seems  that  the  effects  of  the 
compliant  wall  on  turbulence  are  localized,  which  can  also  be  demonstrated  by  comparing 
the  turbulence  statistics  in  figure  7. 

The  magnitude  of  turbulent  kinetic  energy  k,  turbulent  kinetic  energy  production  P,  RMS 
streamwise  vorticity  u'  and  Reynolds  shear  stress  —  (uv)  are  reduced  near  the  compliant  wall 
as  seen  in  figure  7.  The  quadrant  analysis  in  figure  8  shows  that  the  reduction  of  Reynolds 
shear  stress  is  due  to  the  decreased  intensity  of  Q2  and  Q4  events.  As  commonly  observed  in 
turbulent  flows  with  skin  friction  reduction  manipulations,  the  compliant  wall  shifts  the  log- 
law  region  in  the  mean  streamwise  velocity  upwards,  away  from  the  wall  (see  figure  9).  The 
mechanisms  of  skin  friction  reduction  appear  to  be  related  to  the  stabilizing  or  weakening 
effects  of  the  compliant  wall  on  the  near-wall  streamwise  vortices,  which  are  regarded  to  play 
a  dominant  role  in  near- wall  turbulent  transport  phenomena  (Robinson  1991). 

Instantaneous  compliant  wall  shapes  and  the  contours  of  wall-normal  velocity  at  the  wall 
at  two  instants  are  shown  in  figure  10.  The  skin  friction  on  the  wall  is  low  at  one  instant 
and  high  at  the  other.  The  contours  of  streamwise  fluctuating  velocity  in  the  x- 2-planes 
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0.001 


y=+1 


Figure  8;  Comparisons  of  quadrant  analysis.  The  average  position  of  the  compliant  wall  is 
at  y  =  —1.  solid  line:  the  reference  case,  dashed  line:  the  compliant  wall  case. 


Figure  9:  Mean  streamwise  velocity  in  the  compliant  wall  case  near  (a)  the  compliant  wall 
and  (b)  the  rigid  wall. 
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(b)t=2420 


Figure  10:  W^ll-normal  displacement  and  contours  of  wall-normal  velocity  at  two  instants 
of  (a)  low  skin  friction  and  (b)  high  skin  friction  on  the  compliant  wall. 

near  the  compliant  wall  are  plotted  at  the  same  two  instants  in  figure  11.  From  figure  10 
and  figure  11.  it  can  be  seen  that  the  compliant  wall  forms  a  longitudinal  ridge  beneath  a 
high-speed  streak  and  a  corresponding  groove  beneath  a  low-speed  streak,  even  though  the 
contours  of  the  instantaneous  wall-normal  velocity  at  the  wall  are  intermittent  and  oval¬ 
shaped.  [Kang  Sz  Choi  (2000)]  made  similar  observations  for  the  case  of  their  active  wall 
control  for  skin  friction  reduction.  At  the  instant  of  low  skin  friction  on  the  compliant  wall, 
the  near-wall  streaky  structures  are  very  obvious.  Their  breakdown  is  observed  at  the  instant 
of  high  skin  friction  at  the  wall. 

Discussion  First  results  of  our  direct  numerical  simulations  of  turbulence  on  a  compliant 
wall  indicate  a  clear  potential  for  drag  reduction.  While  it  is  true  that  the  results  described 
above  were  obtained  for  a  somewhat  non-generic  case  at  a  very  low  Reynolds  number,  we  find 
that  the  effect  of  the  compliant  wall  is  limited  to  the  near-wall  region,  which  can  be  assumed 
to  show  a  universal  behavior  that  is  only  moderately  affected  by  Reynolds  number.  It  is 
therefore  not  unreasonable  to  assume  that  the  behavior  that  we  found  for  our  low-Reynolds 
number  case  will  persist  into  the  range  of  higher  Reynolds  numbers,  and  for  wide  channels. 
Of  course,  this  assumption  will  have  to  be  tested  in  numerical  simulations  as  well.  We  have 
already  started  such  an  investigation  of  turbulence  in  a  wide  compliant  channel  at  the  time 
of  this  report. 

3.3.2  Conclusions 

The  low-dimensional  model  described  above  indicates  that  it  might  be  possible  to  signifi¬ 
cantly  reduce  the  drag  of  wall-bounded  turbulent  flow  using  a  compliant  wall  of  appropriate 
stiffness.  The  required  stiffness  for  optimal  reduction  of  bursting  frequency  has  to  be  such 
that  the  eigenfrequency  of  the  wall  is  one  to  two  orders  of  magnitude  higher  than  the  bursting 
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(•)t>2360 


Figure  11:  Contours  of  fluctuation  streamwise  velocity  at  the  plane  —  2  parallel  to  the 
average  position  of  the  compliant  wall  at  two  instants  of  (a)  high  skin  friction  and  (b)  low 
skin  friction  on  the  compliant  wall. 


frequency,  which  corresponds  to  a  wall  that  will  show  only  minute  deformations  under  the  ac¬ 
tion  of  the  turbulent  pressure  fluctuations.  This  finding  is  compatible  with  the  experimental 
results  reported  by  Choi  et  al.  [3]. 

Because  of  the  smallness  of  the  wall  deformations  it  is  clear  that  the  internal  damping  of 
the  wall  material  does  not  play  a  significant  role  in  reducing  turbulent  fluctuations.  Rather, 
the  reduction  in  bursting  frequency  seems  to  be  an  indirect  effect  of  the  presence  of  the 
moving  wall,  which  keeps  the  system  close  to  its  fixed  point  for  a  longer  time. 

We  should  also  emphasize  that  we  are  well  aware  of  the  inherent  danger  of  the  simplifica¬ 
tions  we  introduced  in  destining  the  model  that  is  described  above.  It  is  for  this  reason  that 
we  are  currently  concentrating  our  efforts  on  developing  a  direct  numerical  simulation  of  tur¬ 
bulent  flow  over  a  compliant  wall,  which  will  be  used  to  verify  the  predictions  of  our  model 
as  well  as  to  refine  the  model  as  necessary.  In  particular,  once  we  have  the  direct  numerical 
simulation,  we  can  use  the  data  to  directly  calculate  eigenfunctions  for  the  compliant-wall 
case,  so  that  we  do  not  have  to  rely  on  the  ad  hoc  Stokes  eigenfunctions  we  were  using  in 
the  present  study.  We  expect  this  simulation  to  give  us  direct  insight  into  the  modifications 
of  the  flow  dynamics  that  are  caused  by  wall  compliance,  and  thus  to  allow  us  to  refine 
our  modeling  assumptions.  As  a  final  remark,  we  note  that  in  this  project  we  view  direct 
numerical  simulation  and  low-dimensional  modeling  as  complementary  techniques.  Thus, 
our  simulations  will  allow  us  to  improve  the  model,  and  the  improved  model  will  give  more 
accurate  predictions  of  optimal  parameters  for  drag  reduction. 
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•  Blossey,  P.  N.  (1999) 

Drag  reduction  in  near-wall  turbulent  flow.  Ph.  D.  Thesis,  Cornell  University. 

•  Blossey,  P.N.;  Lumley,  J.L.  (1999) 

Control  of  intermittency  in  near-wall  turbulent  flow.  To  appear  in  Intermittency  in 
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•  Rempfer,  D.  (1999) 

Dynamical  Structures  in  Open  Fluid  Systems.  In:  Evolution  of  Spontaneous  Structures 
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•  P.  X.  Blossey  and  J.  L.  Lumley.  Active  control  of  near-wall  turbulent  flow.  United 
Technologies  Research  Center,  East  Hartford,  CT,  January  1999.  Invited. 

•  P.  N.  Blossey  and  J.  L.  Lumley.  Active  control  of  near-wall  turbulent  flow.  Applied 
Mathematics  Seminar,  University  of  Wisconsin  (Madison),  February  1999.  Invited. 

•  P.  N.  Blossey  and  J.  L.  Lumley.  Active  control  of  near-wall  turbulent  flow.  Dynamic 
Systems  and  Controls  Seminar,  University  of  California,  San  Diego,  March  1999.  In¬ 
vited. 

•  J.L.  Lumley  and  P.  N.  Blossey.  Control  of  intermittency  in  near- wall  turbulent  flow, 
Isaac  Newton  Mathematics  Institute,  Cambridge,  UK,  June  1999.  Invited. 
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•  D.  Rempfer  and  J.  L.  Lumley.  Low-Dimensional  Models  for  Near-Wall  Turbulent  Flow: 
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Conference,  Maui,  HI,  June  2000.  Invited. 

•  D.  Rempfer  and  J.  L.  Lumley.  Low-Dimensional  Models  for  Turbulent  Shear-Flow 
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Analysis  of  Turbulent  Flows.  53d  Annual  Meeting  of  The  American  Physical  Society  s 
Division  of  Fluid  Dynamics,  Washington,  DC,  November  2000.  Contributed. 

•  D.  Rempfer.  The  Fable  of  the  Wall-Normal  Pressure-Gradient  Boundary  Condition  - 
Or:  Why  Muenchhausen  was  a  Liar.  53d  Annual  Meeting  of  The  American  Physical 
Society's  Division  of  Fluid  Dynamics,  Washington,  DC,  November  2000.  Contributed. 

•  J.  Gibson.  D.  Rempfer,  J.  Lumley.  Boundary  Conditions  for  Low-Dimensional  Models 
of  the  Turbulent  Boundary  Layer.  53d  Annual  Meeting  of  The  American  Physical 
Society’s  Division  of  Fluid  Dynamics.  Washington,  DC,  November  2000.  Contributed. 

•  K.  Bhaganagar.  D.  Rempfer,  J.  Lumley.  DNS  of  Spatial  Transition  to  Turbulence  in  a 
Boundary  Layer.  53d  Annual  Meeting  of  The  American  Physical  Society's  Division  of 
Fluid  Dynamics.  W^ashington.  DC,  November  2000.  Contributed. 

•  L.  Parsons.  D.  Rempfer,  J.  Lumley.  Low-Dimensional  Model  of  the  Interaction  of  Near- 
Wall  Turbulence  with  a  Compliant  Boundary.  53d  Annual  Meeting  of  The  American 
Physical  Society’s  Division  of  Fluid  Dynamics.  W^ashington.  DC,  November  2000.  Con¬ 
tributed. 

•  S.  Xu.  D.  Rempfer,  J.  Lumley.  Direct  Numerical  Simulation  of  the  Interaction  of 
Near- Wall  Turbulence  with  a  Compliant  Surface.  53d  Annual  Meeting  of  The  Ameri¬ 
can  Physical  Society’s  Division  of  Fluid  Dynamics.  W^ashington.  DC,  November  2000. 
Contributed. 

•  V.  Juttijudata,  D.  Rempfer,  J.  Lumley.  Development  of  a  Localized  Low-Dimensional 
Approach  to  Turbulence  Simulation.  53d  Annual  Meeting  of  The  American  Physical 
Society  s  Division  of  Fluid  Dynamics,  Washington,  DC,  November  2000.  Contributed. 

6.2  Technology  Transitions  or  Transfer 

None. 


6.3  New  Discoveries,  Inventions,  or  Patent  Disclosures 

None. 
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